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^ ; We consider the scattering of quasi-particles off the boundary created during a first 

order electroweak phase transition. Spatial coherence is lost due to the quasi-quark 
damping rate, and we show that reflection on the boundary is suppressed, even at tree- 
level. Simply on CP considerations, we argue against electroweak baryogenesis in the 
Standard Model via the charge transport mechanism. A CP asymmetry is produced 
in the reflection properties of quarks and antiquarks hitting the phase boundary. An 



effect is present at order a^y in rate and a regular GIM behaviour is found, which 
can be expressed in terms of two unitarity triangles. A crucial role is played by the 
damping rate of quasi-particles in a hot plasma, which is a relevant scale together 
with M\y and the temperature. The effect is many orders of magnitude below what 
observation requires. 
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1 Introduction 



A very concise summary of our main ideas and results on Standard Model (SM) baryoge- 
nesis in the presence of a first order phase transition, was recently presented in ref. [|1]. An 
academic^ scenario at zero temperature is developed with all the particulars in ref. 0. The 
aim of the present work is to consider in detail the finite temperature (T) scenario, where 
decoherence effects become essential. 

Gamma-ray and cosmic data do not show any evidence for primary antiparticles on scales 
up to the level of clusters of galaxies. Nucleosynthesis constraints require a baryon number 
to entropy ratio in the observed part of the universe ~ (4 — 6)10~ n at T > 1 GeV 0. 
In the absence of a sensible mechanism for matter-antimatter separation on such large scales, 
a plausible alternative is to require that the microscopic laws of physics are responsible for 
such a baryonic excess [[J. This scenario is called baryogenesis. 

The SM electroweak phase transition occurs rather late in the cosmological evolution, at 
a time when the expansion of the universe is slow compared to weak interaction time scales 
and T ~ 100 GeV. Sakharov's condition^] of departure from thermal equilibrium could be 
fulfilled through a sudden, first order, phase transition f| |§ . In this case, bubbles of true 
ground state, characterised by a non-zero vacuum expectation value of the Higgs field v, 
grow and fill the preexisting v — universe. The evident non-equilibrium element is given 
by the propagation of the bubble surfaces, and physics in their vicinity should be examined 
in view of baryogenesis 0. 

It has been known as well for many years that weak interactions violate baryon number || . 
At low temperatures, the rate is negligible, while at high temperatures the rate is enhanced, 

T u sph ~ n(a w T)\ v = (1.1) 

T s V h ~ exp 



T v^O (1.2) 

where v denotes the vacuum expectation value of the Higgs field, and k is a parameter of 
0(1) or smaller. 

The SM complies as well with the C and CP violation requirements. In particular the 
presence of complex Yukawa couplings induces a phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) fpLOfl matrix for three generations, responsible for CP violation. 

Sakharov's conditions imply that both baryon number-violating and C and CP-violating 
processes have to undergo an out-of-equilibrium period. The currently proposed scenario for 
SM baryogenesis ||11|| is a charge transport mechanism ||12|| . It is graphically summarised in 



Fig. |I|, where we have zoomed into the vicinity of one of the bubbles. In the wall rest frame 
a net baryonic flux is hitting the surface from the unbroken phase. Consider for definiteness 
the trajectory of an initial right-handed quark. Upon eventual reflection, a left-handed quark 
will travel backwards to the unbroken (y = 0) phase. Electroweak interactions have been a 
priori active throughout the trajectory; their CP- violating effects induce a different reflection 
probability for quarks and antiquarks and, thus, a CP asymmetry in the reflected flux. Then, 



4 Academic because the physical first order phase transition is a thermal effect. In ref. we considered 
an hypothetical T = world with two phases of spontaneous symmetry breaking separated by a thin wall, 
to sharpen our tools and disentangle the effects specific to the presence of a wall from the pure thermal ones. 
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Figure 1: Artistic view of the charge transport mechanism, as described in the text. The 
hungry "pacman" represents rapid sphalerons processes. The wiggly lines stand for collisions 
with thermal gluons. Only electroweak loops are depicted, represented by dotted lines. 



unsuppressed sphaleron processes, eq. ( |1.2|) , restoring equilibrium in the unbroken phase, 
can "swallow" the outgoing quarks, transforming the CP asymmetry into a baryonic one. 
The sweeping of the expanding bubble will automatically transfer the latter to the broken 
phase (u ^ 0), where we live. An important survival requirement for the produced baryon 
asymmetry was given in refs. []T3 and ||14j| , by requiring that sphaleron processes inside the 



bubble are weak enough so as not to wash it out. In perturbation theory, this results in an 
upper bound of the Higgs mass ~ 45 GeV, in conflict with experiment. This is a big problem 
for the scenario, although the perturbative treatment may be inadequate, and the question 
is the subject of much work at present. 

We will not enter the discussion on whether a first order phase transition did take place. 
It will be assumed that it did, and that an optimal sphaleron rate is present as well. Our 
aim is to argue, on a quantitative estimation of the electroweak C and CP effects exclusively, 
that the current SM scenario is unable to explain the above mentioned baryon number to 
entropy ratio^l Notice that once a CP- asymmetry on the reflected baryonic current, Acp, is 
obtained, the induced baryon asymmetry is at most n^/s ~ 10 _2 Acp, in a very optimistic 
estimation of the non-CP ingredients P^] fllTf . Acp > 10~ 8 is thus required. 



The symmetries of the problem are analyzed in detail for a generic bubble in ref. || . Major 
unknowns of the scenario are the wall velocity, v wa u, and the wall thickness, I. Typical values 
for \v wa ii\ are non-relativistic, ~ 0.1— 0.4[14| [17]]. Our analytical results correspond to the thin 



wall scenario. The latter provides an adequate physical description for typical momentum 
of the incoming particles \p\ <C l/l- For higher momenta the cutoff effects would show up, 
and indeed "realistic" walls have I ~ 10/T — 20/T, but it is reasonable to believe that the 
thin wall approximation produces an upper bound for the CP asymmetry. We work in a 
simplified scenario with just one spatial direction, perpendicular to the wall surface: phase 



J A discussion of intuitive expectations is given in the introduction of ref. | 
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space effects in the 3 + 1 dimension case would further suppress the effect. It is assumed as 
well that both phases are in thermal equilibrium : this is only broken by the wall motion. 

At T ^ the correct incoming asymptotic states are quasiparticles instead of particles, 
built up of a resummation of the thermal self-energies of the particles. Nevertheless, the 
three building blocks of the problem are essentially the same as in the T = case|2|], 
to wit: the CP- violating couplings of the CKM matrix, the presence of CP-even phases 
associated to complex reflection coefficients for certain values of the energy of the incoming 
quasi-particles, and the well-known fact that, at T ^ 0, the fermionic on-shell self-energy 
cannot be completely renormalized away and induces physical transitions already at the 
one-loop level. Even starting from massless particles, the eigenstates of the effective thermal 
Dirac equation have an "effective mass", or plasma frequency, due to the QCD thermal 
self-energies. It gives the overall energy scale of the problem ~ g s T ~ 50 GeV, for T ~ 100 
GeV. This results in a shift in the position of the reflection coefficient threshold for a given 
quasi-particle. Nevertheless, the tree-level reflection region for a given flavour is still of the 
order of the corresponding current (non-thermal) mass m. In other words, the time required 
for complete reflection of a given quasi-particle is of the order of 1/m. 

On top of the above, the authors of ref.[|TTJ have pointed out that electroweak thermal 

loops are present as well, and because the latter are flavour dependent, their resummation 

can lift the QCD degeneracy of the spectrum of quasi-quarks even far from the wall in the 

unbroken phase! . The result would be further shiftings and reshufflings in the reflection 
p 

coefficient thresholds!!] and, more important, a sizable CP asymmetry due to interference 
between different flavours and chiralities. In our opinion, the important point to retain 
is that CP violation is a quantum phenomenon, and can only be observed when quantum 
coherence is preserved over time scales larger than or equal to the electroweak time scales 
needed for CP violation. 

We argue |J that the quantum phase of the quasi-quarks is in fact lost much before the 
time scales mentioned in the above paragraphs. A further fundamental difference with the 
T = case is the damping rate, 7, of quasi-particles in a plasma. Small momenta are relevant 
for the problem under study, and it is known that at zero momentum, the QCD contribution 
to the damping rate is 7 ~ 0.15 g 2 s T, i.e. ~ 19 GeV. Due to incoherent thermal scattering 
with the medium, the quasi-particle has a finite life-time ~ I/27, turning eventually into 
a new state, out of phase with the initial one. Although the life-time is larger than the 
overall QCD time scale mentioned above, ~ l/(50GeV), it is small in comparison with both 
the electroweak thermal self-energy time scales for any quasi-quark, and with the tree-level 
reflection times ~ 1/m for all quasi-quarks but the top. 

The problem can be rephrased in terms of spatial decoherence. Reflection is a spatial 
property. The quasi-quarks have a group velocity of ~ 1/3 and thus a mean free path, or 
coherence length, of ~ I/67. Over larger scales, such as those needed for instance for total 
reflection, the quantum coherence of light quasi-quarks is damped and any pure quantum 
effect, such as CP violation, is suppressed. 

We thus show that tree-level reflection is suppressed for any light flavour by a factor 

6 Unlike the T = case, Lorentz invariance is lost even far from (without) the wall. At T — 0, far from 
the wall in the unbroken phase, all particles are massless, at any order in the self-energy corrections. Any 
base is a good one, as rotations are physically irrelevant. 

7 We will see that, finally, it just amounts to technical complications. 
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~ m/27. The presently discussed CP- violation observable results from the convolution of 
this reflection effect with electroweak loops in which the three generations must interfere 
coherently in order to produce a CP-violation observable. It follows that further factors 
of this type appear in the final result, which is many orders of magnitude below what 
observation requires and has an "a la Jarlskog" |26| type of GIM cancellations. We show 
as well that the effect is present at order in a perturbative expansion, as in the T = 
case0. 

The scope of the present paper goes beyond the particular issue of baryon number gen- 
eration in the SM. Spatial loss of quantum coherence may be relevant in other microscopic 
processes at finite temperature. A more rigorous treatment than the heuristic one developed 
here would be welcome, though. In our opinion, this general topic deserves further attention. 

In Sect. we discuss on general grounds the relation between baryonic and CP asymme- 
tries. Sect. |3] describes the quasi-particle spectrum far from the wall both in the unbroken 
and the broken phase, as well as at the interface. It contains the dispersion relations, ex- 
emplified by the corresponding effective Lagrangians/Hamiltonians which take into account 
the QCD damping rate. Sect. f| considers the tree-level scattering of a quasi-quark hitting 
the bubble wall from the unbroken phase in a one-flavour world. It parametrises how the 
loss of quantum coherence results in a damping of the reflection properties of the plasma. A 
density matrix formalism is developed as well. Sect. [5| considers the same scattering problem 
in the realistic case of three generations with mixing. It includes the computation of the CP 
asymmetry as well as comparison with literature and comments on wall thickness effects. 
The conclusions are summarised in Sect. || 



2 Baryon Asymmetry Through Charge Transport 

In the introduction, we have discussed qualitatively how it is possible to generate a baryon 
asymmetry in the context of a first order phase transition by fermion transport. In this 
section, we state more precisely the relation between the baryon asymmetry, n^/s, and the 
CP-odd asymmetry that we aim to compute. 

In the fermion transport mechanism, we can distinguish two steps. Firstly, a spatial 
chiral baryon number separation is produced due to the CP-violating scattering of quarks 
on the wall. 

Denote by r™' b , t^ b (f£ b , £!t' 6 ), the reflection and transmission amplitudes from the un- 
broken (u) or broken (b) phases, for an incoming quark (antiquark) of chirality x ^- m our 
convention, a positive velocity is defined as flowing from the unbroken to the broken phase. 
The spin a z is a conserved quantum number. For positive (negative) group velocity, its 
eigenvalue in the unbroken phase is given by chirality +(— )x, as will be shown in section 
( P-2| ). In the broken phase, although chirality is no longer a good quantum number, we keep 
the same notation for the a z eigenvalues. The conserved current density in the ^-direction is 

j z oc u(p')j z u(p) = V z u\p')u{p), for p = p' (2.1) 

8 For anti-quarks, the CP conjugate of a quark of chirality \ has chirality — x- More precisely, the doublet 
quark (antiquark) has chirality —1 (+1), and conversely for singlets. To simplify notations, x wu l indifferently 
stand for (L,R) or the associated 75 eigenvalue (-1,+1). 
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where V is the group velocity. For convenience we will extend eq. (|2.1|) to p 7^ p' } in which 
case it may be taken as a definition of a "non diagonal" V z . For an incoming spinor u x , 
normalised as u"'w" = 1, we have 



x x 



^ = \r u x \ 2 \VyV£\ = \r u x \\ ^- = \t u x \ 2 |Vj/V^|, (2.2) 

Jinc Jinc 

where we used \V™ X /V x \ = 1, consistent with our future approximations. The one particle 
thermal density matrix in the wall rest frame, p x ' b , is a positive definite hermitian matrix i. 
It is not translationally invariant (i.e. it is non-diagonal in momentum), due to the presence 
of the wall and depends on the plasma velocity, —v wa u. Also, it is different for doublet and 
singlet quarks, as electroweak interactions generate an 0(0;^,) difference in their thermal 
self-energies. 

The density matrix of incoming particles for one given chirality \ is normalised as usual, 

dpTi[p u x (p,p)} = l, (2.3) 

where Tr stands for the trace in flavour space. This corresponds to having one particle in 
the total volume L, i.e. to a particle density far from the wall of 1/L. 

The chiral baryon number current density in the unbroken phase, for a unit incoming 
density, is then given by the sum of the corresponding reflected and transmitted chiral current 
densities^, 

Ax = ^ / / dpdp'Tr[P- x (P,P') {r U - x ( P ')r u - x ( P y ~ r>')r» f }] 



+ 



\ V x\l \/l W Wl Tr[^(p,p') {t h x {p')t\ip) S ~ ^ x (p')^x(p)^- ( 2 - 4 ) 



This expression can be simplified by using the relations imposed by unitarity and CPT 
symmetry, or more exactly CP'T. The latter is defined as the combination of CPT and a 
7r rotation around the y axis, which leaves the wall invariant. 

• CP'T 

r x=^x)^ r x ={rl)\ q=(? x y (2.5) 
where superscript t stands for transposition in flavour space. 

• Unitarity, i.e. current conservation, 

«t + \Vi/V-\ t«t«t = !, r ^t + |^ /V J| ^6 1 = L (2 _ 6) 



9 The density matrix for antiquarks of chirality — x is the same, if we assume there are no primeval 
asymmetries. 



In deriving 2.4 we have used the assumption (trivially verified for non-interacting particles prepared 

in equilibrium in a half space) that p! x (p,p') r v L x (p')r u _ x (pf+ p x (p,p')t x (p')t b ' x {p) ] / \V b { P )V b {p')\ 

is the matrix density for the particles of chirality x going away from the wall in the unbroken phase, and 
analogously for antiquarks. The full density matrix further has chirality non diagonal pieces, but these do 
not contribute to the computed current. 
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Then, 

A x = ^ // d P dp'Tr[(pl x - P h x ){rl x {rlj - f^J}] (2.7) 

where the possible flavour non-diagonal elements in {pCf have been neglected, and the p,p f 
dependence omitted for simplicity. 

The total outgoing baryon number current density in the unbroken phase, per unit in- 
coming quark/ ant iquark flux, is: 

B u = A" + Al x = A / jdpdp* Tr[(p« x -pl + p\ - 4){ryryt - f«(f«)t}] (2.8) 

For p u _l x = p x b , clearly B u = 0. This shows that, up to order 0(a w ), only chiral baryon 
number gets separated by the wall. We will neglect this subleading effect and assume p u \ = 

Pjt = P U,b f rom now on - 

Furthermore, neglecting the current (non-thermal) mass effects in the thermal distribu- 
tions 0, p u and p b only differ because of the opposite boost. If p u ' b were the boosted Fermi 
thermal distribution, and for small wall velocities, 

p u -p b ~ v wall fp. (2.9) 

where p corresponds to the Fermi thermal distribution in the plasma rest frame, p oc rip, and 
/ = p~ 1 dp/dv wal1 < 1 0. For the problem under study, p differs from the Fermi distribution, 
but the relation (|2~9| ), with / < 1, should still be a good order of magnitude estimate. 
In this approximation, we define the CP asymmetry as 

A CP = Asign( X ) J J dpdp'Tr[ p K x (r^) f - f^f^}]. (2.10) 

In a second step, the sphaleron transitions are taken into account. We do not consider 
strong sphalerons ||15| , which would further reduce the effect Hl6| , as they require to take 



into account the O(aw) difference in the thermal distributions for doublets and singlets 
mentioned above. Weak sphalerons only affect left chirality particles (or right antiparticles). 
As explained in the introduction, these sphaleron processes have a very different rate in each 
phase, eqs. (|1 . 1| ) and ( p..2|) . In an ideal picture, all the left baryon number in the unbroken 
phase is completely diluted, while the one in broken phase remains. Eventually, when the 
phase transition is completed, the remaining total baryon number, divided by the original 
number n of quarks with one given chirality and spin, is 

nT mal /n = = - l -Al ~ V"A CP . (2.11) 

Only the contributions from the unbroken phase appear, as sphalerons transitions are absent 
from the broken phase. r^ ttmal j s is smaller, as the total entropy s is at least one order of 
magnitude larger than rJli 



n This is a good approximation for quark masses smaller than g s T. 

12 For the Fermi distribution, / = (1 — np) p ^ p , where p u ~ p, p' is the momentum in the unbroken phase, 
and p b in the broken one. 

13 There are 4 spin-chiralities per quark and per antiquark, three generations, and additional quanta in the 
plasma (leptons and gauge bosons). Besides, the quasiparticles correspond to a limited part of phase space, 
E ~ g s T. 
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Unfortunately, the sphaleron processes are not so efficient and computing the actual 
number of sphaleron transitions is a difficult diffusion problem . The authors of give 
an already optimistic estimate for n^/s which is ~ 10~ 2 Acp- Without relying on this 
calculation, it is necessarily true that 

-<A CP . (2.12) 

s 

A C p gives a highly conservative upper bound for n b /s. This will be enough to rule out 
SM baryogenesis within this mechanism, and we will devote the rest of the paper to the 
computation of this asymmetry. 



3 The spectrum of quasi-particles. 

In this section the spectrum of quasiparticles is derived from the thermal loop contribution 
to the quark self-energy. The solutions in a world with just one phase, either spontaneously 
broken (v ^ 0) or unbroken (v = 0), are computed first. Secondly, the case when the two 
phases coexist separated by a boundary (thin wall) is discussed, as well as the matching of 
the two spectra. 



3.1 The quark self-energy. 

This is a rather technical subsection, whose results are needed later. We present a novel 
calculation of the real part of the self-energy, including the electroweak contributions, in the 
broken phase. The computations are performed in the real time formalism. We compare 
to previous results for the unbroken phase. The imaginary part computed in QCD at zero 



momentumflSl is considered. 



3.1.1 Real part in the broken phase. 



We use the notations of generalised to several flavours: 



Re(E(fc)) = -aft -hi- cm (3.1) 

where a, b, c and m are Lorentz invariant matrices in flavour space, u is the four-velocity of 
the plasma and k = (u, k) is the external momentum, m denotes the mass matrix for the 
external flavours. In the plasma rest frame and the mass basis, 

Re(£(a;, fc))7o — —h(u>, k) — a(u, k)a ■ k — c(u, k)m^Q (3.2) 

where 

h(uj , k) = a(u , k) uo + b(u , k) . (3-3) 
The one loop calculation gives the following matrix elements: 

a(u, k) fi = f A{m h 0)6 fi +^[J2 fw,i A{M h M w )+{f z A(m h M z )+f H A{m h M H ))5 fl ] (3.4) 

1 i 
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where the index I runs over the internal flavours for W exchange, and M denotes the corre- 
sponding mass. In this expression, 



f=[h 2 + Q 2 g 2 s 2 w }(L + R) , f m = {(l + ^)L + ^R}K H K* 



if 



If 4 
2^ 



{Tf-Q^f + ^L + i^i-Q^f + ^R , f H = *(L + R), (3.5) 



where L, R are chiral projectors. K represents the Cabibbo-Kobayashi-Maskawa matrix and 
Aj = rrii/M w are related to the usual Yukawa couplings by the following relation, 



fi = 



The integral is given by: 



A(M F ,M B ) = ±f °°^([ 



1_ r oo dp (\ (u> 2 +k 2 +A) p 



Lt( P )-^Lj(p) + ^]n B (E B ) + 



2k E B k J ^ 1 E B i 



where 



and 



Lf{p) = [log 



2kp + 2E B u + u 2 -k 2 + A 
-2kp + 2E B u + uo 2 - k 2 + A 



±[E L 



-Ei 



Lfjip) = -[log 



2kp - 2E F u + uo 2 - k 2 - A 
-2kp - 2E F uo + u) 2 - k 2 - A 



=F [E F — > — E'f], 



with 



E f ,b = ^P 2 + M 2 FjB , A = M| - M|, n FjB (E) = (exp £/T ± l)" 1 . 
The function /i(u;, fc) is given by 



(3.6) 



(3.7) 

(3.8) 

(3.9) 
(3.10) 



h(u, k) fi = -f H(m i: 0)6 fi - Mi H(Mi, M w ) + (/ z #(m,, M z ) + / H H( mi , M H ))6 fi ], 

(3.11) 

where 



The chirality breaking term c(u;, fc) has the following expression 



(3.12) 



c(u,k)fi = 2/ C(mi,0)8fi + 



9 A 



E— ^C(M / ,M W ) + 



(<?z C(m,, M z ) - / H C{ mi , M H ))5 fi 



(3.13) 
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with 



s 2 

and 



9w,i= ^(XfL + X^KuKtf, g z =\[-ZQ i s ^{Tz-Q i s 2 w ) + %{L + R) (3.14) 



C(M F , M B ) = \ f °° ^(£-L+(p)n B (E B ) + £L+(p)n F (E F )). (3.15) 
3.1.2 Real part in the unbroken phase. 



In the unbroken world the expressions (|3.4j )-( 3~T5| ) apply with all masses equal to zero. For 



the sake of comparison with previous literature, we give here the leading (0(T 2 )) contribution 
when u, \k\ <C T. In this limit, the one loop fermionic self energy in the unbroken phase 
has been computed in QCD PH| , |2"I] and generalised to electroweak interactions in JTT]]. We 



find: 

r\ T 2 „ / to \ { 2na s 3iiaw ( tan 2 6 W 1 . , , . . ■-> . , „ 
T 2 / w \ r 2vra s vra^ / 2 2 1 



= -~ F {w\) is + ( <??ta ' 12 ""' + rO} ■ (3 ' 17) 

o(w, fc) L = -h(u, k) L (3.18) 

1*1 b \\k\) 

a(«, k) R = -h(u, k) R ) (3-19) 
cll,R) hi y R correspond to the coefficients of the projectors L, R in a,h, and 

F(x) = Ulog(^)). (3.20) 



with 



and 



2 V \x- 1 



In eqs. (p.!6|) -( 3.19|l it is understood that all terms which do not contain the CKM matrix 
are flavour diagonal^. 

3.1.3 Imaginary part of the Self-energy. 

It is known that the one loop calculation of the imaginary part of the self energy (proportional 
to the damping rate) is incomplete. A leading order QCD computation of the quasi-quark 
damping ratefl9j has been performed at zero k: 7 ~ 0.15g 2 T i.e. ~ 19 GeV at T = 100 
GeV. We will neglect its electroweak component, and assume that in the vicinity of k = 
the damping rate remains close to the above-mentioned value, i.e. 

7 Im(E(w,jfc)) ^ -2 7 (3.21) 



14 The careful reader may notice several differences in numerical factors between eqs (3.16)-( |3.19[ ) and eqs. 
(6.6)-(6.8) in O], once the change of basis is taken into account. 
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3.2 The spectrum. 

The spectrum of quasi-particles in a quark-gluon plasma with only QCD interactions has 
been extensively studied in the literature PD| , [F2TJ. The effects we need for generating a 
baryon asymmetry are electroweak, though, and it is necessary to include the electroweak 
thermal loops. The latter will modify the spectrum in a complicated way, although the 
corrections are quantitatively small with respect to the QCD contributions. This will allow 
a convenient perturbative treatment, as it will be shown. 

Let us summarise first the results in pure QCD, i.e. setting a w = in the self-energies 
( p.2[) . The spectrum of quasi-particles is given by the isolated poles of the propagator (or 



the zeros of the self-energy). At k = 0, it is given by the solution of 

u + hQc D (u,0)- 2*7 = (3.22) 

To gain some insight into the energy scale of the spectrum, we simplify (|3.22| ) neglecting 
all the energy scales other than T (i.e. m ; << T). We work to first order in j/(g s T), an 
approximation we will stick to in the rest of the paper. The solution then is 



Wqcd-T^^ 1 -H = ^qcd-H- (3-23) 

Eq. fl3.23j ) shows that QCD has shifted dramatically the position of the singularities 
of the quark propagator (by singularity we mean the pole and the location of the thermal 
5{p 2 — m 2 )'s), and sets the energy scale of the quasi-particles at 0(g s T) ~ 50 GeV. 

The dispersion relations w(k)qcd satisfied by these quasi-particles are given by the van- 
ishing eigenvalues of the matrix 

ryoS -1 =u + }iqcd(u, k) - 7 5 (a • k)(l + a QC D(u, k)) - 2ry. (3.24) 

We are interested only in the low momentum region, and upon linearising in momentum 
we find, 

w(k) QCD = uj° qcd + x(° ■ k)V QCD ~ iTqcD, (3.25) 

where x is the eigenvalue of 75 (the chirality), i.e. +1(-1) for right (left)-handed fermions, 
and 

_ l + a Q CD(u QCD,0) r 2 7 rq9fi s 

VQCD = 7 77T — , 1 QCD - 7 TvT' {6.2b) 

with 

(, ,0 1 _l d/*QC£>(^,0) | , , 

s Q cd(uj qcd , 0) = 1 + — | u o od . (3.27) 

The functions sqcd are the inverse of the residues at the poles. In a second quantisation 
treatment of the quasi-particles, they can be absorbed through wave function renormalisa- 
tion, in order to properly normalise the kinetic terms. When deriving eq. (|3.25|) higher 



orders in •y/{g s T) have been consistently neglected, leading to real arguments in eqs. (|3.26| ) 
and (ET27D . 



The spectrum is obviously flavour and chiral degenerate, and that equation ( |3.25| ) holds 
separately for each flavour. For simplicity let us take k along the z axis. Vqcd has a simple 
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meaning: up to a sign it is equal to du(k)qcD/dk z , the group velocity of the wave packet. 
From eqs. (|3.16| )- (|3.19| ) it is easy to see that, at leading T 2 order, sqcd = 2 and the group 
velocity is ±1/3. More precisely, for a given value of k z and a given flavour, there are four 
solutions: two possible chiralities and two possible helicities. The group velocity is given 
by the eigenvalue of ~ l/3750" 2 = l/3-y5<j z (k z ) 2 i.e. l/3xhk z , where h is twice the helicity, 
X the chirality and k z = k z /\k z \. As a result, the group velocity has the same sign as 
the momentum when the chirality equals helicity. The solutions with opposite chirality and 
helicity are disregarded at zero T since they correspond to negative energy states. At finite T, 
due to the shift of the energy, part of these branches become physical. They are often called 
"abnormal" branches. As a consequence, for a given value of k z , a normal and an abnormal 
branch have opposite slopes and there is generally one level crossing in the spectrum. It is 
also obvious from the above that positive (negative) group velocity corresponds to 75 = o z 
(75 = — a z ). The two normal (abnormal) branches are degenerate. 

For large \k\ the dispersion relations become non linear and the residues of the abnormal 
branches are exponentially small ~ exp(—\k\ 2 / Uqpp) , \k\ ^> ujqcd, implying that abnormal 
excitations become irrelevant at large momentumEJ. However, we will restrict to the region 
where the linear approximation is appropriate, i.e. for momentum \k\ <C ^qcd- 

In what follows, we consider the electroweak contributions that lift the chiral and flavour 
degeneracy of the spectrum, together with subleading effects. We discuss separately the 
asymptotic spectrum in the unbroken and broken phases, which differ in the quark mass 
effects. In the real situation, with both phases separated by a thin wall interface, it is 
necessary to connect both asymptotic spectra and further approximations will be performed. 
In the following subsections, we start by considering the toy example of one generation in 
order to disentangle the mixing effects and, finally, we formulate the problem for several 
generations. Weak effects are not taken into account in the damping rate, which will then 
be the same for both chiralities and every flavour: 

Tl = T/j = Tqcd- (3.28) 

3.2.1 Asymptotic Spectrum for one generation. 
• Unbroken Phase 

The procedure is the same as in the case of QCD, although taking the functions h(u, k) 
and a(u, k) as given by eqs. ( |3.11|) , (|3.4|) with vanishing masses and a w 7^ 0, for one 



generation. They will be different for right and left chiralities: the degeneracy of the normal 
(abnormal) branches is lifted. We start with the one-loop self-energies in the mass basis, 



V ui + hl + 5- k(l + at) - 2if I x 



15 This, as well as the doubling of the number of poles corresponding to particles carrying a given charge, 
is quite natural when recalling that annihilation operators can act non-trivially on a thermal state, to create 
a physical "hole" . Clearly, for momenta so large that the corresponding state is not thermally occupied, we 
expect a strong suppression. 
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where the index u refers to the unbroken phase, and the functional dependence of h's and a's 
is kept implicit. The on-shell states correspond to the zeros of the determinant of i~fo(S~ 1 ) u , 
and the corresponding eigenstates verify the effective Dirac equation 

ijoiS- 1 )^ = 0. (3.30) 

At k = 0, the solutions are given by the equation 

< >R + h u L , R «,R, 0) - 2*7 = 0. (3.31) 

As we are interested in the low momentum quasi-particles, we expand the effective Dirac 
operator around these poles and linearise in momentum. We also neglect 0(j/g s T), 

HoiS-'T^-^ = [1 + ^jLjJCw - < >R ) - 75(1 + a\ R )a • k, (3.32) 

where we have used the momentum independence of 7. 

The solution for the dispersion relation at small momentum is then, 

w(k) L)R = u u LjR + X (v • k)Vt R - i^ln (3.33) 

where, as before, ui^ R = Re(w^ R ) and 

1 + 0^0(^0,0) 27 

yu = LR\ L,R> J T u 7 

S L,r[ UJ L,R^ U J s L,R\ u L,Ri U J 

with 

dhy R (u,0), 

*5U«*> 0) = 1 + Qijj \»U- (3 - 35) 



The spectrum ( |3.33[ ) is illustrated in fig. 0(a). Even without quark masses, the degeneracy 



between L, R is now broken by the weak interactions due to thermal effects. 

In the regime of low momentum ( |3.33| ), the dispersion relation is well described by the 



following effective Lagrangian for free quasi-particles: 

C eff = ^ L \idt - %V£d -a-u u L )^ L + ¥ R {id t + V%id ■ a - uj\)^ r 

+m L Tl^ L + i^ R r R ^ R , (3.36) 

where ^ l and ^ R are respectively the left-handed and right-handed fields that create quasi- 
particles with the given chirality. £ e // is not hermitian due to the damping rate. 

• Broken Phase 

In this case, the self-energy is given by 

irv ( Q-i\b _( uj + h b R -a-k(l + a R )- 2i 7 (1 - c)m \ . . 

^ b >-{ m t(l_ c )t u + h b L + a-k(l + a b L )-2ij J' {6 - 67) 
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Figure 2: Dispersion relations for quasi-particles in the (a) unbroken and (b) broken phases. 
The full (dashed) lines are normal (abnormal) branches. The upper (lower) lines correspond 
to left (right) chirality. When only QCD loops are considered, upper and lower branches are 
degenerate and intersect the vertical axis at uJq CD . The vertical lines in (b) represent the 
gaps of width ~m, in which total reflection occurs. 

The dispersion relations w(k) defining quasi-particles are given by the zeros of the deter- 
minant of (|3.37|) . In the real part, the effect of the mass is to replace the two lines of the 
unbroken spectrum by a hyperbola. This is again quite similar to the T = scenario, except 
that the center of the hyperbolas is not necessarily at the origin. Between the minimum 
and maximum of the hyperbolas branches, an energy gap of width ~ m appears around the 
center. In this energy gap no plasma excitation can exist. This new spectrum is illustrated 
in fig. (b). 

These very narrow energy gaps are physically crucial since they induce total reflection of 
a quasiparticle coming from the unbroken phase with its energy precisely ranging in this gap. 
However, as we shall argue later, the large value of the damping rate prevents the existence 
of quasiparticles with such a small energy band. 

In order to find the asymptotes we first consider the case k = with vanishing non- 



diagonal terms in eq. Q3.371 ) . The centers of the hyperbolas are given by 

+ 0)=0 (3.38) 

They are slightly different to those in the unbroken phase, real part of eq. ( |3.31|) , due to 
mass effects in h b . The solution with m = 0, \k\ small, is similar to the one in the preceding 



subsection (|3.33| ), with different values of Vl,r and sl,r- A numerical estimate for down 
quarks gives a group velocity V£ = V£ = 0.339 and s b L R (uj b LR ,0) = 1.89, close to the values 
1/3 and 2, respectively, obtained with the unbroken loop at order T 2 . For u and c quarks 
we obtain — V R — 0.346 and h\ R (ou b L R , 0) = 1.88, while for the top the results are 0.165 
and 2.5, respectively. 

In the same way as we did for the unbroken phase, it is possible to obtain the effective 
Lagrangian for quasi-particles in the broken phase for the low momentum regime, 
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C eff = ¥ L {id t - iV b d ■ a - + ^ R (idt + V b id ■ a - u b R )* R 

+M\r L V L + i* R r R V R - {^\ii$>r + Vr^Vl) (3.39) 



where V£ R , T b LR are defined in terms of h b , a b analogously to eqs. Q3.34|) and (|3.35|) , and 



m(l-c^O)) 

4 1/2 K,o)4 1/2 K,o)' 

• The Interface Unbroken- Broken 

The real problem we want to solve is the scattering of quasi-particles on the interface 
between the broken and unbroken phases. The usual procedure in scattering theory in- 
volves matching the in-states (asymptotic states in the unbroken phase) to the out-states 
(asymptotic states in the broken phase) through a unitary scattering matrix. The problem 
is, however, that the spectrum of single quasi-particle states do not form a complete basis, or 
in other words, the residues of all the single quasi-particle states do not sum up to on (Shis IS 
usually associated with collisionless damping which is neglected here altogether. In massless 
QCD (as in NR plasmas), this type of Landau damping only becomes relevant for large 
enough momenta and we expect the same to be true here. . In fact, we have seen that these 
residues s u £ b R are different for L or R quasi-particles and depending on whether they are in 
the broken or unbroken phase. This is a small correction coming from the subleading effects 
in T, so in order to simplify the problem we will neglect them when obtaining the asymptotic 
spectrum. Thus, we will only include in the determination of the quasi-particle dispersion 
relations the following corrections: 

• QCD corrections at leading T. 

• Diagonal Weak Interactions corrections leading in T. 

• Linear terms in momentum. 

This approximation ensures that the residues of the dispersion relations (or wave function 
renormalisation constants) are the same at both sides of the wall and equal for both chiral- 
ities. 

In this approximation, the poles in the unbroken phase are given by the equation 

w\,r + h L , R (w° Li R, 0) - 2*7 = 0, (3.41) 



where hz,R(w, k) contains only the leading-T QCD and diagonal weak corrections of ( 3.11 ). 
Now, expanding r^S^ 1 )" around these poles and linearising in momentum, 

hoiSo'T ^ = [1 + fjl<J(" - ™°l,r) - 75(1 + al tR )a ■ k, (3.42) 
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with aL^(uj, k) containing the leading-T, QCD and diagonal weak corrections in the functions 
( p.4[) . Using eqs. fl3.16|) -( p.l9|) , together with the expansion F(x) ~ 1 + 3^2,2; > 1 we get, to 
leading order in 7/(<7 s T), 



3' 



(3.43) 



1 + 



dh L)R 



2. 



(3.44) 



Notice that although h L R is flavour dependent, the pole condition in eq. ( |3.41| ) implies 
that the above results, eqs. (|3.43|) and ( |3.44|) are not. The damping rate is 
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1 + 



7- 



(3.45) 



The dispersion relations describing the quasi-particles in the unbroken phase are given by 
the zeros of eq. (|3.42|) : 



(3.46) 



In the broken phase, the only difference with respect to eq.( |3.46|) , in this approximation, 
results from the L, R mixing effects due to the tree level mass. Thus, using eq. (|3.42|), 




* 7 o(S -Y = i^ 1 )* + m t * , (3.47) 



where we have neglected the functions c(u, k) because they are subleading in T. The disper- 
sion relations in the broken phase are just given by the determinant of i r yo(SQ 1 ) b . 

From ( |3.42D and ( 3.47Q , we can easily write the effective unperturbed Hamiltonian for 



quasi-particles in the presence of a wall, in the basis of asymptotic states in the unbroken 
phase, 



kcr z d z -i'j + uj 







H° eff =[ ~^~7^ 1 • ^ UKZ \ )• (3-48) 



3.2.2 Several generations. 

In the case of several generations, the exact solution of the spectrum is very complicated, 
as the functions (|3.4j ), (p.ll|) , ( ^.13|) are non-diagonal 3x3 matrices in flavour space. In 
particular, h(uj, k) and 1 + a(u), k) may not be diagonalisable in the same basis, in such a 
way that single quasi-particle asymptotic states cannot be chosen to be diagonal in flavour. 
In order to simplify the problem, we will derive the dispersion relations for asymptotic quasi- 
particle states neglecting flavour mixing and, as we did for the one-flavour case, neglecting 
also the subleading effects in T, both in QCD and weak corrections. As previously stated, 
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in this way the residues of the dispersion relations are the same at both sides of the wall, for 
any flavour and chirality. 

We start with the 6x6 one-loop effective Dirac operator in the mass basis: 



m ' \ u + ht + a-k(l + al)-2i 1 J v ; 



ini t q-v\b _( u + h b R -a-k(l + a b R )- 2i 7 (1 - c)m \ . , 



t(l-c)t to + h b L + a ■ k(l + a b L ) - 2ij 



The results of the previous section hold here, if we substitute all ui° L R and h, a and c by 
3x3 matrices. The expressions for h b , a b and c b (h u , a u and c u ) were given in subsection 
( P-1.2|) with masses different from (equal to) zero. The effective unperturbed Hamiltonian 



for free quasi-particles is then 

H o _ ( -li<Tzd t -ij + u>° R f6(z) \ 

e// " ^ f 0(z) \ia z 8 z - ry + J , ^ 5i J 

where Cj° r and u>£ are 3x3 diagonal matrices (lu° r , to° R , cu° R ) , (uj q1 l , lu° 2 l , cu° 3 l ), respectively. 

In section (^), we will see that it is possible and consistent to include the remaining effects 
(non-leading and mixing) as a perturbation, without modifying the asymptotic states. 

3.3 The damping rate. 

Due to incoherent thermal scattering with the medium, the energy and momentum of the 
quasi-quarks in a plasma are not sharply defined, but spread like a resonance of width 7 [|BJ . 



The QCD damping rate at zero momentum is of the order 7 ~ 0.15 g^T |T]|, i.e. ~ 19 
GeV at T = 100 GeV. Compared to the energy of the quasiparticle, of 0(g s T), the damping 
rate is formally small, i.e. it contains one additional power of g s . This hierarchy allows to 
speak of the mere existence of quasiparticles as coherent excitations 0. 

The quasi-particle has thus a finite life-time ~ l/2 7 , turning eventually into a new state, 
out of phase with the initial one. As the group velocity is of order 1/3, the mean free path 
is of order I/67 ~ 1/120 GeV -1 . Even if it is true that the relative phase between different 
quark flavours is not destroyed by collisions with the QCD thermal bath, it should be stressed 
that the quantum spatial coherence, i.e. the phase relation between points separated by a 
distance > I/67, is lost. 

The imaginary part of the QCD self-energy is much larger than the real part of the elec- 
troweak self-energy. For all quarks but the top, it is also much larger than the corresponding 
mass gap in the broken phase discussed in subsection ( 3.2.1 ). In other words, as we shall see 



explicitly later, the coherent electroweak processes relevant in the following require a time 
much larger than the mean free time of the quasiparticles. 



17 In practice however g s is not so small, and the width of ~ 20 GeV is not really small compared to 
the energy ~ 50 GeV. We could say that quasi particles hardly exist. The plasma is mainly an incoherent 
mixture of states. 
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4 Reflection on the bubble wall: one flavour case. 



In this section we discuss a world with just one flavour, and the tree-level reflection properties 
when a quasi-particle hits the boundary (thin wall approximation) separating the two phases 
of spontaneous symmetry breaking. No electroweak effect is considered, other than the mass 
of the quasi-particle. We quantitatively derive the dramatic effect of the QCD damping 
rate on the reflected density of quasi-quarks. After some considerations in terms of damped 
waves, the system is described in terms of wave-packets which exemplify the mean free path 
and mean free time of the states. We discuss as well the associated density matrix formalism. 



4.1 Waves and wave packets 

Consider the effective hamiltonian, eq.( ft.48| ). As stated above, electroweak loops are irrel- 
evant for the purpose of this section, in which case we could settle uj\ = uj° r = uJq CD (see 



eq.( 3.23 )). The L, R dichotomy will be relevant in the next section. The effective Dirac 



equation describing the quasi-particle interaction with the wall is: 

* - ( -* to V> +i * 0D *a^U»> (41) 

There are many solutions to this differential equation. The problem to solve is the scattering 
of a quasi-quark hitting the wall from the unbroken phase. The initial conditions correspond 
then to a state localised on z < at its creation time, travelling towards the z > region, 
and with a mean free time and mean free path as described in subsection |3.3| . 

No single wave solution fulfills all these conditions. The gist of the problem is reflection, 
though, and waves damped in space and travelling towards the wall should be an appropriate 
heuristic treatment of the initial state, that we proceed to develop here. These solutions have 
then a mean free path (coherence length), but are eternal, that is, have no mean free time. 
All the initial conditions can be accounted for in term of wave packets instead, and we will 
see later that this physically more correct approach leads to the same conclusions. 

An exponentially decaying wave, created in the unbroken phase at a point Zq < within 
a mean free path (I/67) away from the wall, and travelling towards z > 0, is given by 



ifa c {z,t;u>) = e*< t 6(z-z ){ 6(-z) 

+ 6(z)e^[u x , x + r x (n x )u^ x }\, (1.2) 



where x denotes chirality and u X)(T is the normalised Dirac spinor for a particle with spin 
a z /2 and chirality \- n x * s ^ ne incoming momentum of a particle with a z = x, i- e -> positive 
group velocity 

tt x = 3(lu - uo° x + ry) = p x + 3z7, 7r_ x = 3(uj - lu°_ x + 27), (4.3) 
where p x = Re(7r x ). tc x is the corresponding transmitted momentum in the broken phase, 



< = 2^x ~ t-x) + 2V (*x + ^-x) 2 - 9 ™ 2 - ( 4 - 4 ) 
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The reflection coefficient r x (uj) has the simple expression: 

m/2 



with§ 



and 



r x (p x ) = r x {uj) = , (4.5) 

p{u) + e i(t> ^ \p{u) 2 - m 2 /4| 



arg(p(cu) - m/2) + axg{p{u) + m/2) g ^ 



1 u° + 

P(") = g [P X + P-xl = w - x 2 x , (4.7) 

where we introduced the momentum- like variable to stress the analogy with the situa- 
tion at T = ||. The analytic continuation of r(p x ) to the complex plane for the variable 
p x (ui) is straightforward from eq. (|4.5| ). 

For quasi-quarks with current quark masses smaller than 27, 



IWx; 



Tfl 

|r x (^ + z 7 )|<— . (4.8) 



It follows that the reflection probability is suppressed for all flavours but the top. For the 
latter, \r(u)\ < 1, as expected. 

4.1.1 Wave packets 

Our main hypothesis results from the uncertainty principle. The finite life-time (and finite 
mean free path) of the quasiparticles induces a finite spreading of their energy (and mo- 
mentum) spectrum of the order of 27 (67). Wave packets allow a parametrisation of such 
a localisation in time and space. We use gaussian wave packets with a unique width 1/d. 
Our final result will turn out to be independent of this particular choice of packets, the only 
requirement being d < I/37. 

Wave packet solutions to eq. Q4.lD can be expressed as superposition of plane waves. 
Assume that, at time to, an incoming quasiparticle wave packet has been created around the 
position zq < in the unbroken phase, with a central momentum p®: 



/ f\ C "hoc 

*L(*,to;*b,to,Pd,d,x) = / dpe-^-^V^-^u^. (4.9) 



. 2tt 3 / 2 

The choice of spinor with s z = \ ensures the positive group velocity, characteristic of an 
incoming wave. Notice that the eigenspinors of the Dirac Hamiltonian with positive group 
velocity are the same all over one branch of the spectrum (see fig. ||] (a)), as long as we stay 
in the unbroken phase and in the linear regime for the spectrum. 



18 To simplify the notations we use the same symbol r x for the reflection coefficient as a function of ui and 



as a function of p x . The analytic behaviour of eq.(4.5) will be discussed later, see subsection 4.1.2 
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From eq. ( |4.1|) it is easy to derive the time evolution of the wave function ( |4.9| ), in the 
unbroken phase: 



*L(M;^o,£o,Po,d,x) 



= 9(t - t ),/-l-e- " Z " Zn) "^ 3(t " tn)1 ^IMhWt-t.Me^M.HMoL (4 . 10 ) 
V an 1 ' 2 

where the relation between the incoming momentum, p x , and w was given in eq. (|4.3|) . 
Eq ( [4.10D exhibits the motion of the center of the packet with velocity 1/3. The Fourier 



transform of eq. (|4.10|) is 

ty inc (p x ,t; z ,t ,p ,d,x) = 

Q(t - t )^/Z^ e - d2( ^- po)2/2 e-^ 2o e- l( "- l7)( *-* o) M XiX . (4.11) 

The wave function is normalised to 1 at t — to, and its norms decays exponentially as ~ 
e -2-y(t-t ) ^ ue dissipation. Let us now turn to a statistical ensemble of particles. Consider, 
in the unbroken phase, an homogeneous flux of such wave packets travelling towards the 
wall. These particles "decay", i.e. they are scattered, at a rate 2j during their flight. This 
decay must be exactly cancelled by creation at the same rate 27, so that the equilibrium 
distribution of particles can remain stationary. To account for this effect, we introduce a 
continuous differential rate of creation N(z , t ,p )dz dt dpQ. The density of incoming quasi- 
particles of chirality and spin x a t some point z < far from the wall and at time t is given 
by 



2 = /-« d Po too dto /+~ dzoM nc (z, t; z , t Q ,p , d)\ 2 N(z , t , Po ) 

= ^f-™dp N(z Q ,t , Po ). (4.12) 

Taking for the rate of quasi-particle creation 

N(z ,t ,p ) = 2 1 n F (u° x +p /3,T), (4.13) 

the equilibrium Fermi density is recovered, 

~= J dp n F (uJ x + p /3,T), (4.14) 

where np{uS) = l/(exp(u/T) + 1). 

4.1.2 Reflection probability. 

The reflected part of the wave packet, travelling backwards in the unbroken phase is 



/ d r+00 

V? ef (z,t; z ,t ,p ,d) = 6{t -t )\l / _ dp x 
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e -rf , (Px-w) a /2 e -iP-x*-<Px*o e -i(«-i7)(t-to) rx (p x ) u _ X)Xj ( 415 ) 

where the uo dependence is the same than for \I/j nc , as energy is conserved in the interactions 
with the wall, and where r x {p x ) has been defined in eq. ( |4.5|) . In eq. ( |4.15|) , x (— x)represents 
the incoming (outgoing) chirality. 

Let us define the reflection probability n r as the ratio of the reflected flux close to the 
wall on the unbroken side (z — 0~), to the incoming flux. Since the incoming and outgoing 
group velocities have the same absolute values, it amounts to the ratio between the reflected 
density at z = 0~ and the incoming density. Notice that the choice (z = 0~) maximises n r , as 
the reflected flux will be diluted when travelling backwards in the unbroken phase. Further 
than a few mean free paths away from the wall, the reflected flux of particles is transported 
by diffusion rather than free streaming. As we will see in the next section, diffusion could 
be qualitatively described by replacing ( |4.13| ) with its local equilibrium generalisation. For 



simplicity, we concentrate now on the reflected flux close to the wall. 

At time t = 0, the reflection probability is obtained integrating the reflected density, 
eq. ( }4.15|) , over all creation times in the past, to < 0, and over all creation locations in the 



unbroken phase, z < 0, 

\ r+OO rO rO 2 

n*0 = 0,t = 0) = ^ dp dz dt N(z , t ,p ) ** e/ (0, 0; z , t ,Po, d, x) ■ 

^ J — oo J— OO J —CO 

(4.16) 

This is consistent with the creation probability eq. ( |4.12j ). N(zo,t ,p ) was given in eq. 



(Km 



The computation of the reflection probability, eq. ( f4.16| ), requires a careful analysis of the 
mathematical properties of the function r x , eq. ( |4.5| ). r x is analytic on the complex plane but 
for a cut extending on the real axis between p(uj) = —m/2 and p(uj) = m/2. Close to the cut, 
~ 7r/2 in the upper half-plane, while ~ —tt/2 in the lower half-plane. When integrating 
on the real axis we will take the determination from the upper half-plane. Outside the cut 
on the real axis and for p{uj) > m/2, = while for p{uj) < —m/2, = 7r0. It results that 
|r x (u;)| < 1. More precisely, |r x (o;)| = 1 for —m/2 < piyj) < m/2, and |r x (o;)| < 1 outside 
that segment. This is illustrated in fig. (|3[). Fig (g) displays the phase of r x (u), and shows 
that it varies very fast over the real axis on the cut. It is worth to notice that this cut does 
correspond to the total reflection domain for plane waves. It is a very small domain for light 
quarks. 

In eq. ( 4. 15| ) the integral over p x is strongly reduced by the rapid rotation of the phases of 



the exponential and the function r x . This is at the origin of the dramatic suppression of CP 
asymmetry that we advocate in this paper. We have checked this suppression numerically, 
although the same result follows from a very simple analytic argument, which we proceed 
to develop. The only singularity in the integral dp x in eq. Q4.15D is the cut on the real 
axis of the function r x . Consider a rectangle in the complex p x plane with sides given by 
: the real axis, a parallel to it in the upper half-plane, and two vertical sides at ±oo. The 
extension of the integral in eq. ( 4.15| ) to this contour vanishes since no singularity is present 



inside the box. The integrals on the vertical sides vanish due to gaussian suppression. The 
integral in eq. ( [4.15|) on the real axis can then be replaced by an integration on a parallel 



19 Then, on the real axis, r x — m/2/ (p(u) + ^p{uj) 2 — m 2 /4) for p(uj) > m/2, and r x — m/2/(p(ui) — 
x/piio) 2 — m 2 /4) for p(ui) < —m/2. 
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Figure 3: Plot of \r(p)\ on the upper complex p(u) half-plane, with m = 0.5. The lower 
left edge is the real axis. The crest at \r(p)\ = 1 on this axis corresponds to total reflection. 
Notice the rapid decrease of \r(p)\ with increasing imaginary part. 



axis, J^+iiJj^ dp x , with < (3 < 1. The condition (3 > ensures that the integral over z in 



( f4.16| ) is finite, while the condition (3 < 1 does the same for the integral over t . 

The choice (3 = 1, and thus Im(p x ) = 3^7 is equivalent to change variables p x 
in eq. (QgjB 



Px + 3i 7 



■$* ef (z,t; z ,t ,p ,d) 



0(t - 1 ) 



2vr 3 / 2 



dp x e 



- rf2 (""x-Po) 2 /2 -iir- x z-iTT x z -iu>(t-to) 



(4.17) 



The reflection coefficient r x (ir x ) is then bounded as displayed in eq. (fO). The exponen- 
tial in (|4J.5|) also provides a factor exp(9(i 2 7 2 /2), and a bound for ^* e f results as long as 
d« l/^ 



^ ef (z,t;z ,t ,p ,d) 



< 



m 



4 Jp(u° + Po/3) 2 + 7 2 + 0(m 2 



0(3d 7 ), 



(4.18) 



where p(u) has been defined in ( |4.7| ) and the bound in ( |4.8| ) has been slightly improved for 
later use. 



20 Notice that p x + 3i7 is the variable tt x defined in eq. (|L3|). 

21 Physically, this inequality means that the layer of particles sitting "on" the wall has a negligible thickness 
d, compared with the layer of particles that ever have a chance to reach the wall before being damped. 
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Figure 4: Plot of the phase ofr(p) on the upper complex p(uj) half-plane, with m = 0.5. The 
phase is zero on the real axis outside the total reflection domain, —0.5 < x < 0.5, and varies 
rapidly inside this domain. This explains the suppression observed when the integration is 
performed on the real axis. 



From eq. (|4.16 ) it follows that 



9 

m 



n r (0,0)< . (4.19) 

167 7 2 + 0(m 2 )S 



It is worth to notice that the results expected for zero damping rate 7 = can be recovered. 
In this case, eq. ( |4.19|) reduces to 

m 

n r (0,0)~—, (4.20) 

This is indeed the result that would follow directly for vanishing damping rate, when the 
reflection probability can be directly computed in terms of plane waves: 

-J dp x n F (u,T)\r x (uj)\ 2 ~- (4.21) 



where u and p x are related as in eq. 

The comparison of eqs (|4.19| ) and fl4.20|) shows a ~ m/ (I67) suppression of the reflection 
probability when the damping rate is taken into account, due to the fast rotation of the phases 
in the integral in eq. ( |4.15|) . This is a typical quantum effect: the different frequencies in the 
integral, eq. ( |4.15 ), add up as waves, not as probabilities. We may express the phenomenon 



at work as follows: the damping rate destroys the quantum coherence which is necessary to 
have a large reflection probability, through an average over a rapidly rotating phase, 
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< r x >< r % ><< |r x | 2 >. (4.22) 

The spread in energy and momentum of the quasi-particles results in a smearing of the 
reflection integrals, contrary to the case of vanishing damping rate, where the states are 
monochromatic, and thus well represented by plane waves. 

4.1.3 A simple approximation to the reflection probability. 

A useful simplification is obtained if the integration on to in eq. ([4.1 6|) is extended to +oo. 
A qualitative estimate of the error introduced by this trick is given here. A rigorous proof 
is developed in appendix [A]. 

The presence of the potentially dangerous exponential e 7 * in the expression for ^ re f, 
eq. Q4.15| ), does not lead to any divergence, due to the gaussian smearing. Indeed, the only 
effect of extending the to integration from to +oo is to add a gaussian tail. The width of 
the gaussian being d, the added tail will be 0(exp(— l/(7<f) 2 ), i.e. small as long as dj <C 1. 
More precisely, as explained in appendix 0, the correction is negligible for 

3^7 < 1. (4.23) 

The reflection probability (|4.16|) is then given by 

27 rO poo d 

n*(0,0) = ^y dpo J dz J dt n F (uJ x + po/3,T)-^-^ 

/oo r 00 

dp x / ^e- d2( ^- po)2/2 - d2( <- po)2/2 e i < 2 °-^ 20 e i( "- w ' ) ' r x (7r x )r*(4). (4.24) 
-00 J—oo ' v 

The integral over t simply gives 2tt5(uj — u>') = 6n5(p x — p' x )- The z dependence is then 
exp(— 67^0), whose integral is trivial. Finally, 

n*(0,0) = - / dp n F (u; x +po/3,T)— / dp x e- d2 ^-^ 2+ ^ 2 \r x (n x )\ 2 . (4.25) 

£1 J—oo y TV J—oo 

This is the main result of this section. After neglecting the to > terms and upon the 
change of variables p x — > p x + 3«7 and correspondingly u — >• u> + ij, the reflected flux for 
wave packets is just a gaussian smear-out of the reflection fluxes for plane waves, as naively 
expected. 

Two comments are appropriate. The first one concerns particle number, which must be 
conserved. Unitarity implies that the damping of reflected flux is exactly compensated by 
the enhancement of the transmitted one. The underlying physics can also be understood. A 
very localised particle penetrates the broken phase in a time I /3d, whereas for small masses 
it requires a rather long time 2/m to be significantly reflected (the reflection dominantly 
happens in an energy window of gap ~ m, the total reflection domain, and by uncertainty 
principle this corresponds to a time ~ 1/m). If in between the particle is "measured" 
in the broken phase by collisions with the plasma, the wave function collapses into one 
corresponding to a localised particle in the broken phase, and as a result, the particle has 
been effectively transmitted. 
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Secondly, the d dependence of ( |4.15|) can be clarified. When the mass is large enough, 



7 <C l/(3d) <C m/2, reflection typically occurs between 2 collisions, and the effect of both 
the 7-shift and the 1 / (3<i)-smearing can be neglected, nicely recovering the usual asymptotic 
reflection probability for plane waves. For the opposite limit of small masses, the d depen- 
dence can also be mild. Indeed, exchanging the order of integration in ( |4.25| ) leaves us with 



a thermal average of the reflection probabilities |r x (7i~ x )| , according to an effective thermal 
distribution 



n d F (u° x + Px/3, T) = 4= f°° %>e- d2(p *- po)2+(M7) + Po/3, T). (4.26) 

\/TT J-oo 



If tif is smooth on scales 1 / (3d) (as is the case, close to equilibrium, if 1/ (3d) <C T) , ~ rip 
and the average ( |4.25| ) is close to a thermal average of monochromatic reflection probabilities. 



We have checked numerically this d independence of our result. If the distribution prob- 
ability was arbitrary and far from equilibrium (i.e. far from np(u)), our linear response 
approach to the relaxation processes loses its meaning, and finer knowledge of all relaxation 
processes becomes necessary. Let us thus stress that we assumed the quasi-particle distri- 
butions to be close to equilibrium (and effectively reducing the reflection coefficient makes 
this assumption easier to meet, for a given wall speed), and moreover required 

7 « « T & g 2 s T « « T, (4.27) 

which is always satisfied for weak coupling g s <^il. Our model for decoherence accounts for 
the physical consequences of a non-vanishing damping rate. In real QCD g s is not small, 
and we believe that a more realistic treatment would lead to a further suppression of the 
CP-asymmetry than advocated in this paper. Indeed, as the damping rate increases the 
description of the system in terms of quasi-particle ( an essential ingredient of the approach) 
looses its physical meaning. 



4.2 Density matrix. 

In this subsection, the results of the preceding one are reproduced using the formalism of 
density matrices. The language of density matrices is very convenient for plasma physics. It 
will clarify how the suppression effect found above depends on the quantum aspects of the 
density matrix, and their relation with the breaking of thermal equilibrium in the process 
under study. 

The effective Hamiltonian in the plasma rest frame, eq. ( |4.1|) , can be split into hermitian 
and non-hermitian parts: 

H eff = H herm - i 7 . (4.28) 

We consider the reduced density matrix p for one quasi-particle at equilibrium in the unbro- 
ken phase. It obeys the following Boltzman equation, 



d t p(t) = -i[H herm , p(t)) +G-L. (4.29) 

G (L) is a gain (loss) term corresponding to the creation (disappearance) of quasiparticles 
due to collisions with the medium. To simplify the problem, we used a linearised equation. In 
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this approximation, the L term is related to the damping rate by L = — 27 p. At equilibrium, 
the probability loss due to L must be compensated by the G term. However, and this is 
our essential assumption, G creates quasi-particles from the medium which are out of phase 
with respect to the ones destroyed by L. We further assume G to be homogeneous in space 
and time in the unbroken phase. A simple example is given by the creation rate defined in 
eqs. ( |4.12|) -( [4.13 ). The resulting equation for the reduced density matrix of a quasi-particles 



with positive group velocity (incoming) and chirality x, is then given by: 

d t p{p,p';t) + i[H herm ,p(p,p';t)} = -2^p(p,p';t) + J dz 6{-z ) 

dp N(z , t, Po)ty in c(Pi t; t, po, zq, d, x)*L(p'> *l Po, zo, d, x) • (4.30) 



S has been defined in eq. ( 4.12Q , L is the volume of the one dimensional box in which the 



density matrix is normalised, and the 9{— zq) factor in the inhomogeneous term accounts 
for quasiparticle creation in the unbroken phase. \I/j nc was defined in eqs. ( (4.9[ )-( OT ). The 
solution of eq. (|4.30|) is 



27 f 

p(p,p';t) = — J dz dtodp 9(-zo)9(t - t )n F (uj + xPo/3,T) 

%inc(P, t] t , PO, Zq, d, X)^inc(p': *l *0, Po, z 0, d, x) , (4.31) 

where eq. Q4.13] ) has been used. From eq. ( |A.7| ) in appendix ^ and eq. ( |4.14| ) it follows that 



dpTT[p u x (p,p)] = l. (4.32) 

It is possible to see that eq. ( [4.31] ) reduces to the equilibrium distribution in the absence 
of the wall, i.e., when the factor 9(—z ) is taken off. Integrating then z from — oo to +00 
yields a 2ir5(p — p') factor, consistent with translation invariance in the absence of a wall, 
and p becomes diagonal in momentum as expected. The to-integral cancels the 27 factor in 
eq. flPlD , 

p{p,p;t) = — J dp n F (u x + po/3iT) • ( 4 - 33 ) 

In the infinite volume limit, 2tc5(0) ~ L, and the po integral gives a gaussian smearing of 
the equilibrium distribution, 

1 r e ~d 2 (p-p ) 2 ~ n F (uj T) 
p(p, P'' t ) = zJ d Po nF ( u x + Po/3, T) ^ 1/2 ^ rfT>1 ^ — (4-34) 

since the smearing can be neglected for d large enough, i.e. 

dE 1 

T»A£ = — A Po ~— , (4.35) 
dpo id 

which is the usual condition ( |4.23| ), ( |4.27| ), and the Fermi equilibrium distribution is obtained. 

In the present problem, the wall at z = requires to treat separately the Zq < and 
z > regions, though. Only the first one is involved in describing the incoming flux, and the 
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Zq integration in eq. Q4.31 ) provides an extra non-diagonal piece arising from the principal 
part of l/(p — p'). This term is exponentially suppressed when \p — p'\ > 1/d, as (|4.11| ) 
imposes that both p and p' lie within a distance 1/d from some p . 1/d is then an upper 
bound for the non-diagonality of p. Depending on the shape and height of the wall, the 
zq > contribution may ultimately reduce the maximum \p — p'\, as it was shown above 
when the wall was removed. This non diagonality of p is a quantum effect as is manifest in 
the oscillation of the exp (— i(p — p')(t — t )/3) factor present in ( |4.31| ), through ( |4.11| ) and 

A simpler way to describe the effect follows by computing the Fourier transform of p 
with respect to (p — p'), in order to obtain the Wigner function f(p, z) p3f . A sharp cutoff 
imposed in z§ then translates into inhomogeneities of f{z,p) on scales at least bigger than 
d, the individual particle diameter. The fact that f(z,p) is not a real positive function (as 
it would be in the semiclassical limit) is a sign of quantum effects. 



4.2.1 Reflection probability from the density matrix. 

The thermal average of any physical observable O is given by Tr(pO)/Tr(p). In particular, 
the reflection probability is the thermal average of T x (p x )r x (p' x ), with r x {p x ) as defined in 
( f4.5| ). For a diagonal density matrix, i.e. in the absence of quantum coherence, the reflection 
probability amounts to the average of \r x (p x )\ 2 , and this is what happens for vanishing 
damping rate (see the discussion in subsection ( [4.1. 2| )). 

The density of incoming particles at position z < is given in terms of the density matrix 

by 

d mc {z) = f f dpdp'e t{p ~ p ' )z p(p,p') = \ if < -~ < 1. (4.36) 
Z7T J J L d 

In fact d mc (z) depends only slightly on z when 37c? <C 1: 

'Ztwall 1 

d inc (0) = =— ~ - if 3 7 rf < 1, (4.37) 



as can be seen from eqs ( |A.7| ) and ( |A.9| ) in appendix 0. For the problem under study, with 



a non diagonal density matrix, the reflection probability is given by the ratio of the density 
of reflected particles over the incoming density, i.e. 

L r 

n r = — I dpdp' p(p,p']t = 0)r(p)r^(p'), (4.38) 



2tt 

where time translational invariance was used to choose t = 0. As expected, eq. ( [4.38[ ) is 
equivalent to ( |4.16|) and the calculation will proceed exactly as in the preceding subsections. 



5 Several flavours: CP violation. 

In this section, we include the non-diagonal weak contributions to the effective Hamiltonian 
in a perturbative expansion. We then compute the CP-asymmetry and compare our results 
with existing literature. 

22 A diagonal density matrix is just a probability distribution, i.e. it is purely classical, it exhibits no 
quantum coherence. 
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5.1 Effective Hamiltonian 



Only the quarks with the same charge are mixed, and the Hamiltonian will be factorised for 
down and up sectors. In section ([3]), we found the effective unperturbed Hamiltonian ( 3.51[ ) 
in the case of several generations, which we recall here: 

It contains only the leading T, QCD and flavour diagonal electroweak corrections. It has 
been linearised in momentum. In order to obtain an observable asymmetry, the remaining 
corrections, i.e. the non-diagonal weak corrections and those subleading in T, are essential 
ingredients. The full effective hamiltonian is 



H eff = H° eff + 9{-z)5H u eff + 6(z)5H b eff , (5.2) 



where 



c H u,b _ 1 ( [6h% (u, k) - i5a^ b (u, k)a z d z ] ^mc(u,k)9(z) \ , . 

eff ~2\ \mc{uo, k)9{z) [Sh^iu, k) + ida u /(LU, k)a z d z ) ) ' 1 j 

Sh u ' b = h u ' b — h and 5a u,b = a u ' b — a are flavour dependent and non-diagonal, as defined 
and discussed in subsections (|3.1.2j) , Q3.2.1|) and (|3.2.2j ). They have to be also linearised in 



momentum. To be consistent with our approximation (i.e. that these effects do not modify 
the residues derived from the free hamiltonian ( |3.48| )) we must expand these functions around 



the poles (|3.31|) and take only the zero order. This means that corrections of 0(k5h, a, c) are 
neglected, which is a reasonable assumption for the low momentum region we are considering. 
Finally, the hamiltonian we must consider is: 



n e ff - u eff + 2 v{ z)y Q Q) j + 2 0{z) y mc ^^ Q) Sh b L ^o t o) J , 

(5.4) 

where we can take u° = uJq CD , i.e., it contains just the dominant QCD part. The pre- 
cise choice of u° in the arguments of the perturbed functions only influences subdominant 
contributions to the asymmetry. 

5.2 Reflection Matrix and CP- Asymmetry 

We are now ready to compute the CP-asymmetry defined in Sec. 2, eq. ( |2.10| ), 

A CP = ^sign( X ) J Jdpdp'lp TrK x V« x - (^)V^}]. (5.5) 

Being a CP-odd asymmetry it must contain the contribution of CP-violating phases, for 
which at least three generations are needed. In this case, the reflection and transmission 
amplitudes r, t are 3x3 matrices in flavour space. 

We show below that an effect is present al order in rate, through subleading T 
corrections. At leading order in T, 0(T 2 ), there is no contribution at this electroweak order, 



28 



because the only flavour dependence is through Yukawa couplings. The physical reason is 
CKM unitarity. The order effect is tantamount to a one-loop electroweak correction 
in the amplitudes. Denote by % and / the external initial and final flavours, and I (or 
I') the internal one in the loop. Convoluting two of these diagrams gives a contribution 
to the rate oc Im{K*fKn(K*, * Kyi)*}. If the function accompanying this factor contains no 
antisymmetric dependence on the internal flavours, the contribution will vanish when the 
sum over I and I' is performed. At leading order in T, the dependence on Yukawa couplings 
is insufficient to produce a non-zero effect, because they factorise as a symmetric function of 
I, A C p would then be a higher order electroweak effect, i.e., a^- Subleading T corrections 
do contain an explicit internal quark mass dependence, which results in the above mentioned 
antisymmetric functions. 

In section 4, we have shown that the reflection probability of quasi-particles in a single 
flavour world can be written as a gaussian smear-out of the reflection probability for planes 
waves (with zero damping rate) after making an analytic continuation p x — > p x + 3ry, eq. 
( |4.25| ). This result holds as well in the case of several generations: 

<(0, 0) = - / d Po n F (u° x + po/3, T) 

" f °° ^ x e- rf2 ^^) 2+ ^ 2 Tr[r x (vr x )rt (vr x )], (5.6) 

71 J-oo 



where now r x (7r x ) are 3x3 matrices and Tr refers to the tracing over the flavour indices. 
The problem is then reduced to calculate the reflection matrices for plane waves r x (p x ) = 
r x (u), evaluate them in the complex point tt x and, finally, perform the gaussian and thermal 
averages. 

Let us start by computing r x (uj). For 5h = 0, the problem reduces to three independent 
one-flavour cases, whose solution we know, eq. ( f4.5|) . An analytic expression in the case 
of the simplified effective Hamiltonian, eq. ( |5.4|) , is still too difficult to find. We will use a 
perturbative approach to obtain an analytic result. The procedure involves deriving an exact 
equation for the reflection amplitude, and proceeding then by iteration, in powers of Sh^^R. 
This corresponds to an expansion of the asymmetry ( |2.1(J| ) in a w sin8 c , as the CP-violating 
effects are contained in the non-diagonal weak terms. 

We work in the flavour basis that diagonalises the unbroken effective Hamiltonian, -ff"//- 
The transformation to this basis from the one corresponding to the unperturbed Hamiltonian, 
H®jp has the form 



I* = M? X. !,x U<¥, (5.7) 





where ip x are spinors with three flavour components, J3 is the identity matrix and Ol is a 
unitary 3x3 matrix. In the prime basis, the solution for right quasiparticle with energy uo 
coming from the unbroken phase has the usual form 



il% e {z,t ] u) = e- iut { 6(-z) e^ z u x , x + e-^ z r x ( Px )u 



-x,x 



+ 6(zy^[u x , x + r x (p x )u„ x J , (5.8) 
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where p x is the incoming particle momentum as defined in eq. (|4.3|) . 

It is possible to obtain an implicit equation for and tl by imposing, 

z > ] H b eff O ip' = u ip'. (5.9) 
For an incoming right-handed quark, eq. (|5.9|) gives 

r' R [u R - u + \sh b R (u°)] + [co° L ~ w + \sh h L {u°)W R 

+ K' R [ 1 - c(uj°) }mr' R + \m[l- c(u°) ]t = 0, (5.10) 

where r' r = OlTr- As Ol is unitary, Aqp is unchanged upon replacing r — > r'. 
In the same way, the equation for antiparticles reads 

f' R [u R l -{8h R {^)Y\ + [Cj° l -u, + \{8h h L {u°)y]f' R 

+ l -f' R { 1 - c(u°r \mf' R + \m[l- c(u Y f = 0. (5.11) 

A similar equation for L reflection amplitudes can be written, although it is not needed to 
compute A C p. 

We look for a perturbative solution of equations ( |5.10| ) and (|5.11|) in powers of a w : 

r R = r ( R ) + r { R ] +r ( R ) + ... (5.12) 

r R is then the solution in the limit Sh,Sa,c = and it reduces to (|4.5| ), as expected. To 
lighten the notation, we will skip the subindices R on r's from now on. The first and second 
orders are given by 

(!) _ 1 ( r^5h b{ R ] + 5/i 6 £V(°) - r^5c {1) mr^ - m5c ' 



2 d{j 



(5.13) 



and 



r 



(2) 



1 {r^5h bi R ] + 6h b( i ] rW _ r ^Sc {1) mr^ - r^5c {1) mr^ 



2 d{j 



1 ( r (i) mr (i) + r (°)Sh b{ R ] + 5h b ? ] r^ - r^5c {2) mr^ - m5c {2) 



(5.14) 



2 dij 

where the indices refer to flavour and dij = {oj 1 l + u R — 2a;) + + mj 2 r " . 

It follows from the discussion in subsection 4.1 that r° are complex in the region of total 
reflection, and this region is the same for particles and antiparticles. This gives rise to CP- 
even phases. In order to generate an observable CP-asymmetry, they must interfere with the 
CKM one, contained in the non-diagonal weak corrections: 5h b R , Sh b L , c. Using eq. (|3.11|) , 
and neglecting those terms in Sh^ that have the same dependence on internal masses (or 
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Yukawa couplings) than Sh R , because they give a zero contribution at this order , we can 
easily obtain: 



Sh b R = a w Xi\ f K u K^I R {Mf) , Sh\ = a w ^ K H K* f I L (Mf) (5.15) 
i i 

and 

c=^J2 K ^MMl), (5.16) 

where we have defined 

I R (M?) = -^H(M h M w ), I L (Mj 2 ) = \ 2 I R (Mj i ), I m (M?)=Tr\iM l C(M h M w ). (5.17) 

It then follows that the first effect in the asymmetry appears at 0(a^) and it comes only 
from the interference of the 0(a w ) effects in Sh R and Sh L . Consequently, there is no effect 
at 0(a 2 ) at leading order in T, because at this order Sh b R = 0. It is interesting to analyze 
the expression for the non-integrated asymmetry at this order, where the GIM mechanism 
is explicitly operative: 



A^p = Tr[ r-W + V*- - 1 + r^V*- 2 -* — antiparticles } 
~ £ Im[ ShlhSh'kU x Im{rr [ + > {{ «> + f{±- + ±- ) ] }. 

(5.18) 

(2) 

A^ P can be shown to have the following structure: 

Agp ~ a 2 (2iJ) T int T ex \ (5.19) 

where J, T mt and T e:ct contain the expected "a la Jarlskog" behaviour of the asymmetry as 
a function of the weak angles (J), the internal quark (T mt ) and the external quark masses 
(T ext ). The connection between Q5.18| ) and ( |5.19| ) is 

Im[Sh b L ) jt 5h b R h} = a&\ i \ j 2i y £lm[K tt iq i K Vj K;, i ](\t - \ 2 V )I R {M 2 )I R {M?) 

i,v 

= a 2 w A i A i (±2iJ)T in *, (5.20) 

with 

J = ±Im[Ki i K? j Ki, j K* H ] = cic 2 c 3 s 2 S2S3S 5 , 

and 

T int = ^ (A 2 _ A 2 +i)/fl(JW? )I R ( Ml+1 ). (5.21) 
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Figure 5: (a) The triangle in the complex plane, with area T mt . (b) The triangle in the 
complex plane, with area T ext . In both cases, the sign is negative (positive) if the arrows turn 
(anti)- clockwise. This figure is only illustrative, a realistic triangle would look extremely flat 
due to the huge mass hierarchy between quark masses. 



The ± sign in J refers to the cyclic order of i,j and I + 1 must be understood as modulo 
3. T mt is the oriented area of the triangle formed by the three points in the complex plane 
given by \fln(M?) + iln(Mf), corresponding to the three flavours of the internal quarks: 
1 — 1,2, 3, shown in fig. (a)) From this construction, it is obvious that the triangle would 
collapse into a line whenever two internal masses are degenerate. 
Finally, 



r 



1 ™i((r&) 2 - (r?.-) a ) r% f l 1 



Although this expression is more involved, it vanishes whenever two external quarks 
masses are degenerate. This can be seen explicitly in a triangle construction for the analytic 
continuation to u + ij, eq. ( |4.25| ), in the limit m << 7. Then, 

TTl 

dij ~ d + O(-) = -2(17 + u) + uj° r + oo° L , (5.23) 

7 

leading to 

rnext Mo I + 2Re(d ) x ^ , o o* 0^,0* \ /r 0A \ 

1 u~TZ x i x i+n r i+ii+i r u - r iiT i+ii+i), (5.24) 

l«o| i 

where again the sum over i is modulo 3. T ext is proportional to the oriented area of a triangle 
with vertices in the complex plane given by A^r^, Fig. ^)-(b). When any two external quarks 
masses are degenerate, two of these points would coincide and the area would vanish. This 
shows that the GIM cancellation also works explicitly for external quark masses, as expected. 



5.3 Results and Comparison with literature 

We now give the numerical results we have obtained at 0(a4) in the reflected baryonic flux, 
both for up external quarks (u,c,t) and downs (d,s,b). We have used the following values 
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20 30 40 50 60 70 80 



Figure 6: The dominant (0(a^)) non-integrated CP asymmetry, Acp(u), when mass effects 
are included inside thermal loops, as a function of the quasi-particle energy, u>. The figure 
corresponds to charge 2/3 flavours. 



for the masses in GeV, Mw = 50, Mz = 57, = 0.006, m s = 0.09, m;, = 3.1, m u = 0.003, 
m c = 1.0, m t = 93.7. The couplings are X d = 1.2 10~ 4 , \ s = 1.8 10" 3 , \ b = 6.2 10~ 2 , 
\ u = 6.2 10~ 5 , A c = 2 10~ 2 and A* = 1.88, and a s = 0.1, a w = 0.035. Fig. (§) shows A(u) 
for the ups (the dominant contribution). The result for the averaged ensemble asymmetry 
is, 

Aqp = 1.6 10~ 21 , Aqp = -3 10" 24 . (5.25) 

In both cases, we find that the asymmetry is dominated by the two heavier external quarks. 

These results are several orders of magnitude in disagreement with the previous results 
obtained by Farrar and Shaposnikov in ref. [IT|. Their estimate is Acp ~ 10 ~ 8 (see their 



eq. (10.3)). This tremendous discrepancy can fortunately be understood. There are mainly 
two important effects these authors did not consider. The first one, which explains the orders 
of magnitude difference, is that they disregard the effect of the damping rate in the scattering 
of the particles on the wall. In our calculation, however, it is an essential ingredient. The 
second difference, responsible for a small factor in the discrepancy, is due to the fact that 
they just considered the leading-T corrections in the self-energies. We have seen in the 
previous section that the subleading terms contain the dependence on the internal masses, 
essential for CP- violation. When they are neglected no 0(a 2 ) effect can appear, and the 
first contribution would loom at O(a^). 

In order to prove that these latter considerations are correct, we have reproduced their 
results, i.e. we have computed the CP-asymmetry just considering leading-T corrections in 
( |5.4|) and setting 7 = 0. The result is shown in Fig. (^-a) and it is in complete agreement 
with the results in [Dj (see Fig.5 in this reference). Improving their numerical calculation 



by the inclusion of the damping rate, gives the result shown in Fig. (0-b). The integrated 
asymmetry from fig (0-b) is ~ 4 • 10~ 22 , to be compared to the result in fig (0-a) of [11], 
~10" 8 . 

With or without the damping rate, it is also possible to work out an analytic expression 
at third order, from the perturbative calculation ( |5.12| ),( |5TT8| ). For 7 = 0, the expected GIM 



cancellation appears, although there is no simple formula that explicitly shows them. On 
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47.9 48 48.1 48.2 48.3 48.4 20 30 40 50 60 70 80 

(a) (b) 

Figure 7: (a) shows the non-integrated CP asymmetry (Acp) produced by down quarks in 
the narrow energy range which dominates for zero damping rate, when masses are neglected 
in the internal loop, (b) shows the dramatic effect of turning on the damping rate effects, in 
the same approximation. 

the other hand, in the case 7 ^ and in the limit m « 7 §, the expression for the peak 
value of the asymmetry beautifully reduces to 



a max 
t-^CP 




J 



(m 2 - m 2 c )(m 2 t - m 2 u )(m 2 c - m 2 u ) (m 2 b - m 2 )(m 2 - m 2 d )(m 2 b - m 2 d ) 



(27)' 



(5.26) 

This was expected from naive order-of-magnitude arguments. 

Finally, the results ( 5.25 ) show that non-leading effects in T give the main contribution 
to the asymmetry in the case of non- vanishing damping rate and, in contrast with [11|, the 
up-sector dominates the asymmetry. 

Very recently, Huet and Sather|^| have analyzed the problem. These authors state that 
they confirm our conclusions. As we had done in ref . [Jj , they stress that the damping rate is 
a source for quantum decoherence, and use as well an effective Dirac equation which takes it 
into account. They discuss a nice physical analogy with the microscopic theory of reflection 
of light. They do not use wave packets to solve the scattering problem, but spatially damped 
waves, as in our heuristic treatment at the beginning of Sect. HI. 



5.4 Wall thickness. 

Notice that the derivation in sect. |] is totally independent of the shape of the function 
r(k). The only requirement was a singularity structure limited to a cut in the region of total 
reflection. This is quite generic: only for very special wall shapes can other singularities be 
expected. For instance, when the wall is not monotonous, a pole with an imaginary part 
may express the decay of a quasi-bound state trapped in a potential well. 

The thin wall approximation used in this paper is valid only for wall thickness I <C I/67, 
while perturbative estimates suggest / > .lGeV -1 > I/67. The CP asymmetry, generated in 

23 This is valid for down external quarks, the case we considered 
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the presence of such more realistic walls, would be orders of magnitude below the thin wall 
estimate, eq. (|5.25| ), reinforcing thus our conclusions, because a quasi-particle would then 
collide and loose coherence long before feeling a wall effect. 



6 Conclusions. 

At finite temperature a plasma is mainly an incoherent mixture of states. The scattering of 
quasi-quarks with thermal gluons induces a large damping rate, 7. Quantum coherence is 
only maintained between two inelastic scattering processes of this type. It is lost over spatial 
distances larger than ~ I/67 ~ (120 GeV) -1 , and over time scales larger than ~ I/27 ~ (40 
GeV)- 1 . 

We have considered both an heuristic approach to quasi-quarks in terms of spatially 
damped waves, and a more physically sound one in terms of wave packets which incorporates 
the above mentioned characteristics of coherent mean free time and coherence length. With 
these tools, we have shown that coherent tree-level reflection of quasi-quarks hitting the 
bubble wall is suppressed for all flavours but the top, by a factor m/2j. This effect might 
also be relevant in certain non-relativistic systems, e.g. in solid state physics. We have shown 
as well that a CP-violation reflection asymmetry appears already at order aw in amplitude, 
and it is suppressed by further powers of m/27. It fails to explain the baryon asymmetry by 
more than 12 orders of magnitude. 

Any non-standard electroweak scenario for baryogenesis, where quantum coherence is 
required over distances and times larger than the above mentioned ones, is subject to the 
same flaw. In particular, "realistic" electroweak walls are thick, much larger than the co- 
herence length. No (quark) coherence should be required through the whole wall thickness 
in a successful scenario. In a more positive view, it is worth to stress that promising candi- 
date theories are those where CP violation is associated to the top and/or heavier flavours. 
Leptonic induced baryogenesis may also be safe. 
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A Proof of the simplified formula for the reflection 
probability. 



This appendix develops a bound to the error introduced in eq. ( 4.24 ) when the integra- 
tion f^clto in eq. ( |4.1_6| ) is replaced by f*™ dt . We derive an upper bound of v)— x,X 
\l/ r e/(0, 0; z , t , p 0j d, A comparison of eqs. ( |4.9|) and (|4.15| ) shows that: 



uL XtX ^ ref (z,t; z ,t ,Po,d,x) 
1 



27T *'* 



dp x e ip -* z r x (p x )y inc (p X} t; z , t , p , d, x) 



As p_ x = p x + 3(u> x — ), eq. ( |A.1[ ) can be expressed as a convolution 



ul^refiz, t; z , t , p , d, x) 



u x>x / dz x e 3i ^x u -x>^ inc ( Zl ,t;z ,t ,po,d,x)r x (z + z x ) 



where 



'2tt 



-ip x z 



The proof goes as follows. We first show that 



-iJ-9 (z) / d9 

V 7T J-tt/2 



cosfm^ sin 9) cos 2 9 



where m is the quark mass. This gives the bound 



(A.r 



(A.2) 
(A.3) 

(A.4) 



\r x {z)\ < \j-m0{z). 



(A.5) 



Next, this bound is inserted in eq. ([A.2Q and the final bound on the to > contribution to 
the reflection probability in eq. ( [4.16[ ) follows. 



Bound on f x (z). 



The integral in ( |A.3[ ) is computed using (|4.5|) and the analytic properties of r x (p x ), 
explained after eq. (|4.16| ). For z < 0, the contour can be closed on the upper half 
plane, with a vanishing result. Else, the cut has to be crosses with the result: 



d Px e-^ z r x ( Px ) = 6{z) / d Px e~^ 



Px+ i Y m2 -Px PX~ i \J m ' 2 -p 



(A.6) 



/ d9 cos(mz sin 9) cos 2 9 

2 J-tt/2 



where a change of variables, p x — > msin9, was performed. This completes the proof 
of eq. ( [A.4| ), and hence of eq. ( |A.5|) . 
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The factor u^— x,X takes care of the spinor u~XtX i n ^re/j so as to keep only the spatial dependence. 
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Integrals on z and t > for the incoming wave packet. 

Before considering the effect of extending the integration on t in eq. ([4.16 ) to +00, 



we discuss the simpler exercise of integrating the squared norm of the incoming wave 
^inc- The particle density 5 in eq. ( 4 . 1 2| ) has been computed far from the wall in 



the unbroken region. For that reason, the integral on z extends to +00. Consider 
the same density close to the wall, i.e. for z = 0~ (and t = to be specific). Since 
only packets created in the unbroken region are relevant, we must now constrain the 
integral on zq from —00 to 0. This is the density of incoming quasi-particles hitting 
the wall, E walL . 



•—•wall 



dp dt dz 6(-z )e(-t )\^f nc (0, 0; zq, t , p , d)\ 2 N(z , t ,p ) 

< 2~ / dp ° N ( z °> to ' = S - ( A ' 7 ) 

where the inequality stems from the positivity of the integrand and the comparison 
with the integral in eq. (|4.12 ). It is useful to notice that, if the integral over to in the 



first line of eq. ([A. 7 ) is extended to +00, we recover 5 i.e. the same result as extending 
the integral over zq to +00. Let us estimate the difference between 'E wal1 and H, i.e. the 
effect of the integral over to > on the incoming density. Afterwards we will extend 
the result to the reflected probability. 

From the last line of eq. ( |4.10| ), it follows that 



S - ~ wal1 = ^ dpo dx Jf dt e-^ +2 ^N(z , t , p ) (A. 
= 1^ dp n F (u° x + p /3, T) dxe-i (e 6 ^ - 1) 



where the change variables to x — —z + 1 /3 was performed, and eq. ( |4.13| ) was used. 
Expanding to first order in 37c?, it follows that, 

H-S tta "~^H«3 if 3 7 rf<l (A.9) 

7T ' 

where eq. ( |4.14| ) was used. The difference between S and "E wal1 consists in gaussian 
tails of size ~ d, small compared to the range of damping I/37. The same phenomenon 
is present in ^ re f, through ^i nc , which appears in the convolution ( |A.2| ). 

Integrals on zq and to > for the reflected wave packet. 
From eqs. ( |A~2|) , flCIOD and (gj, it follows 



\^ref(0,0;z ,t ,p ,d, X )\<-^-^J o d Zl e 53= (A.10) 
The contribution to the reflection probability, for to > 0, is then bounded by 
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2 

m 



n*(0, 0; t > 0) < — / eZ^^K + p /3) x 

/ dzidz 2 dz dt e ^ +27 *° (A. 11) 

Jo 

where a change of variable z — > — z was performed. With the supplementary variables 
changes 

Z\ + z + t /3 = ydcos9, z 2 + z + t /3 = ydsiia.9, z + t /3 = xd, t = t , (A. 12) 
the integral over t leads to 

[ 3Xd dt e 2lt0 = — (e^ xd - l) ~ 3xd } (A.13) 
Jo 27 v ' 

while the integral on x, assuming sin 6 < cos#H, leads to 

/ dx3xd 2 = —{ysind) 2 (A. 14) 

Jo 2 



and, finally, 



It follows that 



3d 4 r dyy 3 e~^ [* d6 sin 2 6 = 3d 4 - — -. (A.15) 
Jo Jo 4 



m 2 d 3 6'j /it — 2\ . . 

n*(z = 0, t = 0; t > 0) < \—) ■ (A.16) 



When this result is compared to eq.( [4.25| ), 



2 

m 



n*(0,0)~ — (A.17) 
7^ 



it is clear that the t > contribution is indeed suppressed, 

{ldf{Ed) < 1 (A.18) 



<(0,0;t >0) . ,. 2 _ 



7^(0,0) 

where we have used the fact that H = 0(g s T) (a one-dimension density of quasiparticles 
with energy in the 0(g s T) range), and that it is easy to choose d such that the condition 
( |1.27p is met simultaneously with (7<f) 2 (S<i) <C 1. For example, 1/d ~ g s T satisfies all these 
criteria. 



25 By symmetry, the alternative sector, sin 9 > cos 8 gives the same contribution, leading to a factor of 2 
in the final result. 
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